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|H ' Small galaxies consisting entirely of population III (pop III) stars may form at high 

C/3 ' redshifts, and could constitute one of the best probes of such stars. Here, we explore 

, the prospects of detecting gravitationally lensed pop III galaxies behind the galaxy 

cluster J0717.5+3745 (J0717) with both the Hubble Space Telescope (HST) and the 

^f) I upcoming James Webb Space Telescope (JWST). By projecting simulated catalogs 

. of pop III galaxies at z w 7-15 through the J0717 magnification maps, we estimate 

' the lensed number counts as a function of flux detection threshold. We find that the 

I ongoing HST survey CLASH, targeting a total of 25 galaxy clusters including J0717, 

|0 i potentially could detect a small number of pop III galaxies if ^ 1% of the baryons in 

^— ^ ' these systems have been converted into pop III stars. Using J exposures of J0717, 

Tij" I this limit can be pushed to ~ 0.1% of the baryons. Ultra-deep /FFS'T observations of 

' unlensed fields are predicted to do somewhat worse, but will be able to probe pop 

' III galaxies with luminosities intermediate between those detectable in HST/CLASH 

, and in JWST observations of J0717. We also explain how current measurements of 

^ ■ the galaxy luminosity function at z = 7-10 can be used to constrain pop III galaxy 

•j^ \ models with very high star formation efficiencies (~ 10% of the baryons converted into 

/\ . pop III stars). 
H ■ 
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1 INTRODUCTION 

The first, metal-free (population III) stars are predicted to 
form in isolation or in small numbers within dark mat- 
ter minihalos of mass 10'^-*' Me at z 10-50 (e.g 
Tegmar k et~ai]|l997l: lYoshida et al.ll2003l : iTrenti fc Stiavellil 
2009; S tacy et al.ll2010h . Some of these population III stars 
(hereafter pop III) likely explode as supernovae and en- 
rich the surrounding medium with heavy elements, initiating 
the transition to the population II/I star formation modes 
known from the low-redshift Universe. 

As more massive dark matter halos (10^~* Mq) start to 
assemble, small galaxies capable of sustaining prolonged star 
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form a tion probably emerg e at 2 « 10-15 (e.g. iGreif et al.l 
l2008l : [Tohnson et al.ll2009f ). Since detailed simulations sug- 
gest that most of these galaxies should form in high-density 
regions that were pre-enriched by pop III stars in minihalos, 
these galaxies are not expected to be completely metal- free 
(e.g. iGreif et all |2009| ) and are most likely dominated by 
chemically enriched stars. True pop III galaxies may, how- 
ever, still form in low-density environments where Lyman- 
Werner feedback has prevented prior star formation in mini- 
halos, and into which sup ernova eject a from adjacent regions 
has n o t yet migrated (e.g.lScannapieco. Schneider fc Ferrara 



liooi; iTornatore, Ferrara fc Schneidej l2007l : iTrenti et al 
[2009'; StiavcUi & Trcnti 2010). 

Simulations of the collapse and fragmentation of metal- 
free gas generically predict that pop III stars should 
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form with high chara cteristic masses (for a review, see 
iBromm fc Larsonll2004l ). but the exact initial mass function 
of these objects remains poorly constrained. About a decade 
ago, many studi es favoured pop III masses on the ord er of 
~ 100 Mq (e.g. iBromm et al.lll999l : lAbel et al.ll2002l ). but 
the most recent simulations hint at masses closer to ~ 10 Mq 
and an initi al mass function that may even extend below 
^ 1 M a (e.g. Stacv ct al. 2010|: lciark et al.ll201ll : lGreif et al] 
l2011al )a: iHpsokawa ct al, 201ll ). Since isolated pop III stars 
are likely to be outsi de the reach of even the James Webb 
Space Tel escope!^ fe.g. [Gardner et aLll2006l : lGreif et al ]|2009l : 
iRvdberg e t al. 201 2) . at least before they go s upernovae 
l|Weinmann fc LilhIBoOSI : IWhalen fc Frroilboiol ). pop III 
galaxies may offer one of the best routes to observationally 
pinning do wn the prope rties of pop III stars. 

While llnoud (|201lh and IZackrisson et all (|2011ah have 
argued that pop III galaxy candidates can be singled out 
from James Webb Space Telescope {JWST) multiband sur- 
veys because of their unusual spectra, the intrinsic lumi- 
nosities of these objects may be too low for detection in 
even the lon gest JWST exposures (| Johnson et al.l l2009l : 
IStiavelli fc Tre nti 2010). Pop III galaxies at z « 6-10 can, in 
princ i ple, also be ident i fied in HST imaging surveys (|lnouel 
I2OIII : IZackrisson et al.1 l2011bh . although this places even 
harsher demands on their luminosities. Here, we will ex- 
plore to what extent gravitational lensing may brighten the 
prospects for detection. 

Lensing by foreground, low-redshift galaxy clusters 
can boost the fiuxes o f high-redshift objec t s by factors 
of u ~ 10-100 (e.g. iBradlev et al.1 [ 2OO8I: IZheng et af 
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2009; 'Brada c et aD 12 009: Mai zv et al.l . 

2009 , ,2011 ; Hall et al.|[2012 ; Bradlev et al.ll2012 ; Zitrin et al.. 
2012D . but at the expense of zooming in on a high-redshift 



volume that is smaller by the same factor. Because of these 
competing effects, surveys targeting lensed fi elds are not 
equa l ly suitable for all sou rce populations (e.g. iMaizv et al.l 
I2OI0I : iBouwens et~al]|2009l 'l . 

In this paper we investigate whether high-redshift {z > 
7) pop HI galaxies behind the 2 = 0.546 galaxy cluster 
MACS J0717.5~f3745 (hereafter J0717) may be detected 
with either the HST or JWST. This cluster has the largest 
angular Eins tein radius detect ed so far (~ 55" for a source 
at Zs ~ 2.5: IZitrin et al.ll2009h . Due to the complex merg- 
ing nature of this cluster, it has a relatively shallow surface 
mass-density profile which boosts the p rojected area over 
which high magnifications are attained (jZitrin et al.ll2009l . 
I2OIII ). Both of these properties combine to make this object 
one of the best lenses currently known (quite possibly the 
best) for the study of faint objects at very high redshifts. 
Previously, we have demonstrated that z ~ 10 dark stars 
could in principle be detectable in a JWST im aging survey 
targeting this object (jZackrisson et al]l2010ah . 

Our simulations are presented in Sect. [2] In Sect. |3l 
we explore the prospects of detecting lensed pop HI galax- 
ies at > 7 in a hypothetical JWST survey of J0717. The 
merits of this approach are also compared to those of ul- 
tradeep JWST observations of unlensed fields. Sect. |4]con- 
tains the correspondin g analysis for the HST/CLASH survey 
l|Postman et al.ll2012l ). which targets a total of 25 galaxy 



clusters, including J0717. Constraints on the properties of 
pop HI galaxies based on existing observations of the galaxy 
luminosity function at z = 7-10 are presented in Sect. 14.21 
A number of uncertainties in our analysis are discussed in 
Sect. [5] 



2 MOCK SAMPLES OF POPULATION III 
GALAXIES 

The galaxy cluster J0717 at z = 0.546 (see Fig. [TJ is a 
very unusual object, with an unrelaxed morphology, a cor- 
respondingly shallow density profile and the largest Einstein 
radius so far detected (Zitrin et al. .2009,). The area over 
which at 2: > 6 background objects attain high magnifica- 
tions /i is exceptionally large: « 3.5 arcmin^ for fi > 10. This 
potentially makes it the best lensing cluster currently known 
for studying faint objects at very high redshifts. J0717 is one 
of the targets of the ongoing 

HST/CLASlE 

survey, covering 

25 galaxy clusters in 16 filters, and also an exciting target 
for upcoming studies with JWST. 

To explore the prospects of searching for pop HI galax- 
ies in a highly magnified field like J0717, and to weight 
the merits of this approach against observations of unlensed 
fields, we generate mock pop III galaxy catalogs for both 
of these situations. The computational machinery for this 
involves the following steps: 

(i) For the lensed samples, simulated pop III halo cata- 
logs (see Sect. 12. ip are projected through the cluster mag- 
nification maps (Sect.[2]2}. This produces the positions and 
magnifications of all potential pop III galaxy images within 
the J0717 field. For the unlensed samples, the pop III halo 
catalogs are simply left as they are. 

(ii) Each distinct halo in these catalogs is assigned a ran- 
dom stellar population age (Sect. [23]). 

(iii) Using the spec tral synthesis model Yggdrasil 
(Zackriss on et all l2011al ), we then derive the HST and 
JWST fiuxes of each pop III galaxy image within these 
mock samples, as a function of various model parameters 
(see Sect. 12. 4) . The fluxes in the lensed samples are adjusted 
by a factor equal to the magnification of each separate im- 
age. Due to the image splitting associated with strong lens- 
ing, a single pop III galaxy can appear at several locations 
within the J0717 field, each with a different magnification. 



2.1 Population III halos 

We construct Populati on III halo ca t alogs based on the nu- 
merical simulations of iTrenti et al.l (|2009l '). To summarize, 
we run a high-resolution cosmological simulation with dark 
matter only {I = 7h^ Mpc, N — 1024'^ particles, a mass reso- 
lution of 3.4 X lO^A^Q and a force resolution of 0.16/i^ kpc), 
which we use to construct halo catalogs and their merging 
histories. We then follow the chemical enrichment of the gas 
within those halos, coupling analytical recipes for cooling 
and star formation (plus subsequent metal production) with 
a wind outflow model that assumes propagation of metals 
out of Population II galaxies with a fiducial velocity of 60 km 
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Figure 1. 4 X 4 arcmin magnification map for tlie largest len s 
currently known, the J0717 galaxy cluster l lZitrin et al.l [20091) . 
The color bar on the right shows the log]^Q(/i) magnification values 
plotted in the figure, assuming a source redshift of z ~ 10. 



s~ . Halos that arc flagged as chemically pristine and meet 
the requirement for cooling via either molecular hydrogen 
(in presence of a self-consistent Lyman- Werner background) 
or Lya, are assumed to be able to form metal-free stars. 
The minimum hal o mass for cooling is shown in Fig. 1 of 
iTrenti et al. Ii200^, and is not repeated here. At 2 < 13.4 it 
corresponds to a dark matter halo with virial temperature 
T^„ = 10* K. 

As discussed in IStiaveUi fc Trentil (|2010l ). the star form- 
ing metal-free halos at z < 15 prefer low-bias regions, and 
form preferentially away from larger Population II galax- 
ies which would otherwise pollute them via wind enrich- 
ment. The predictions, for the formation rate of metal- free 
star-forming halos from our model are shown in Fig. 4 of 
iTrenti et all ([.2009'), and ag ree with similar studies, for ex- 
ample those carried out bv iTornatore. Ferrara fc Schneid"eil 
l|2007l ) if one assumes in our model a relatively low efficiency 
of formation of Population III stars, for a total of a few 
10^ Mq to a few 10^ Mq per halo (containing ^ 10^ Mq in 
baryons). However, given the uncertainties in the processes 
that lead to the formation of metal-free stars, we here ex- 
plore a much wider range of star formation efficiencies, as 
discussed in Sect. 12.41 

In Fig. [5] we plot the number of pop III galaxies per 
arcmin^ and unit redshift throughout the range z = 7-20, 
assuming that such galaxies lose their pop Ills identities 
after Tpopiii — 10^ yr (see Sect. l2.4l for a more thorough de- 
scription of this parameter). The surface number densities 
of pop III galaxies peak at 2 ~ 10, but remain high (> 100 
arcmin"^) throughout the interval z ~ 7-14, which is suf- 
ficient to expect lensing surve ys to be able to catch a few 
(as previously pointed out bv IStiavelli fc Trentil I2OI0I ). To 
see why, consider the following example. The image-plane 
area for which 2 « 10 sources are expected to attain magni- 
fications yU = 10-30 (with an average magnification Jl ~ 17) 
is ~ 2 arcmin^ in the J0717 field shown in Fig. [U whereas 
the corresponding area for fi = 100-300 (/I ~ 160) is ~ 0.2 
arcmin^. A rough estimate of the number of 2 ~ 10 pop III 
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Figure 2. The predicted number of pop III galaxies per arcmin^ 
and unit redshift in an unlensed field, assuming that these galaxies 
lose their pop Ills identities after Tpopiii = 10^ yr (see Sect. 12.41 
for a more thorough description of this parameter). 

galaxies with ^ — 10-30 in the J0717 field would then yield 
300 objects arcmin"^ x 2 arcmin^ x (1/17) ~ 35, whereas the 
corresponding number of galaxies with /i = 100-300 would 
give 300 arcmin"^ x 0.2 arcmin^ x (1/160) ^ 0.38. Hence, 
the surface number densities of pop III galaxies should be 
sufficiently high to allow detection of objects with magnifi- 
cations up to /i < 100 in this field. 

The total masses of these pop III galaxy halos are in the 
range ~ 10^-10* Mq, with likely baryon masses an order of 
magnitude lower. Since JiyST surveys of unlensed fields will 
be unable to detect pop HI galaxies wi th stellar population 
masses below ~ 10^ Mq (|Zackrisson et al. 2011a), a fraction 
> 0.01-0.1 of these baryons must be converted into pop 
HI stars to push them above the detection threshold. In a 
survey of a lensed field behind J0717, the required fraction 
can potentially be lower. 

2.2 Gravitational lensing 

In the process of producing mock samples of lensed pop III 
galax;ies, simulated pop III halo catalogs (see Sect. 12.11) are 
projected through J0717 magnification maps. For this pur- 
p ose, we use the mas s (and magnification) model published 
m IZitrin et al1(|2009l ). which can be consulted for further de- 
tails. The adopted image-plane resolution of these maps is 
chosen identical to the HST /ACS resolution, with 0.05 arc- 
sec pixel" ^. Since this cluster shows uniquely high magnifica- 
tions, and over a large area, we work with a 100 times better 
resolution (0.005 arcsec pixel"^) in the source plane, so that 
magnifications of up to = 100 ca n be credibly tracked . 
While this paper was in preparation, iLimousin" et al.l (|201ll ) 
presented new spectroscopic data on a number of lensed 
background objects in the J0717 field. Even though we have 
not revised our model or calculations to take advantage of 
these new data, their likely impact is deemed to be minor 
(see Sect.l5.4|). 

To facilitate the lensing calculations, the pop III halo 
catalogs are projected onto source planes with thickness 
Az — 1.0, centered on 2 = 7, 8, 9... 15. While the simu- 
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lations predict a few pop III halos at even higher redshifts, 
such objects are out of bounds of the surveys considered, and 
are therefore not treated any further. The simulation boxes 
are comparable to the HST and JWST fields of view in the 
plane of the sky, but much too small in the line of sight di- 
rection to fill up the light cone volume between two adjacent 
source planes. Several boxes are therefore randomly rotated 
and stacked to produce each source plane catalog. This pro- 
cedure admittedly scrambles some of large-scale clustering 
of pop III galaxies, but since these objects display very little 
clustering ( Stiavelli fc Trentill2010l ). any bias resulting from 
this is likely to be small. If anything, the final results should 
be conservative, since the probability of having rare, un- 
usually massive pop III halos in a given lens plane becomes 
smaller when simulation cubes are recycled, compared to the 
case where the entire line of sight volume is self-consistently 
simulated. 

For practical reasons, our lensing scheme necessarily in- 
volves a number of simplifications. Magnifications fi > 100 
are not considered, since the resolution limits of our lens 
and source maps makes the exact magnification difficult to 
assess for such cases. While our lensed object maps typi- 
cally include a few images of objects at z ^ 7-15 in this 
magnification range, all such images are excluded from fur- 
ther analysis. Moreover, all pop III galaxies are treated as 
point sources. A more realistic treatment could imply that 
different parts of the sources exhibit slightly different mag- 
nifications. Since most of the light from pop III galaxies may 
stem from nebular emission rat h er than direct star ligh t 
(e.g. ISchaereil I2OO2I : llnoud I2OIII : IZackrisson et all l2011al l. 
and since the nebulae are likely to extend over a larger area 
than the star-forming regions in these objects, differential 
magnification could alter the effective ratio of nebular to 
stellar emission in sources seen through J0717. The exact 
effect of this is difficult to assess at the current time - the 
spatial extent of pop III galaxies likely depends on the ages, 
stellar initial mass function, halo masses, star formation ef- 
ficiencies and Lyman continuum escape fractions of these 
objects, and simulations covering the full parameter space 
are not yet available. However, unless these pop III galax- 
ies are l arger than the lensed z ~ 7 galaxies detect ed so 
far fe.g. IZheng et al.ll2009l : fBradlev et al1l2008l . \201^ . dif- 
ferential magnification is unlikely to significantly distort the 
observables. 



2.3 Monte Carlo simulation of Pop III galaxy ages 

Within each lens plane, ages are randomly assigned to ev- 
ery halo using a fiat age distribution betw een t = and 
21 Myr (the highest age considered in the lTrenti et al.ll2009l 
simulations). These ages represent the time since the onset 
of pop III star formation within each halo, and are used in 
conjunction with spectral synthesis models (see Sect. 12. 4p 
to derive the fiuxes of the pop III galaxies. For every com- 
bination of model parameters (see below), this procedure is 
repeated over and over to generate a large number of mock 
pop III samples for statistical analysis. 



2.4 Spectral synthesis 

To generate the intrinsic spectral energy distributions 
(SEDs) of the pop III galaxies, we use the Yggdrasil spec- 



tral sy nthes is model ([Zackrisson et al.|[2011al ) with lSchaereij 
12003 ') and iRaiter et al.l (|2010l ') stellar population spectra 



for pop III stars. The nebular contribution to the overall 
SEP is computed u sing the photoionization code Cloudy 
(jFerland et al.lll998h . assuming a spherical geometry for the 
photoionized gas, a filling factor of 0.01 and a constant hy- 
drogen density of n(II) = 100 cm"''. Various model grids 
(SEDs and fluxes in various filter systems) results are pub- 
licly available from the lead author's homepag4f|. 

When creating mock pop III samples, a number of dif- 
ferent model parameters are considered. These are: 

• T he pop III stell a r initi al mass function (IMF). 

As in IZackrisson et al] (|2011al V three di fferent sc e narios 
are considered: pop III.l, pop III. 2 and a iKroupal (|200ll ) 
IMF. The pop III.l option corresponds to a very top-heavy 
IMF with a power-law IMF (diV/dM a M~") of slope 
Q = 2.35 throughout the mass range 50-500Mq . Pop III. 2 
corresponds to a slightly less extreme IMF, with log-normal 
shape, characteristic mass Mc = 10 Mq, dispersion o" = 1.0 
and w ings extending from 1-500 Mq . Finally, the iKroupal 
(j200ll) option results in an IMF identical to that observed 
for more metal-rich stars in the local Universe. 

• Lyman continuum escape fraction, fasc- A sigiufi- 
cant fraction of the ionizing radiation produced by the pop 
III s tars may escape dire ctly into the intergalactic medium 
(e.g. I Johnson et al.ll2009l). with potentially strong effects on 
the resulting SED (IZackrisson et all 1201 lal ). Here, we con- 
sider the two limiting cases /esc = and /esc = 1. The 
/esc ~ scenario result s in a SED dominated by nebu- 
lar emission (type A in IZackrisson et al.l l2011ah . whereas 
/esc ~ 1 (type C) results in a purely stellar SED. Because of 
their radically different SEDs, the colour signatures of these 
two types of pop HI galaxies differ significantly, as discussed 
in Sect. 15.51 Since fasc ~ implies a significantly higher lu- 
minosity per stellar mass, such pop III galaxies are much 
easier to detect. 

• Star formation history. Two options are considered, 
either an instantaneous burst of star formation (giving a 
single-age population) or a constant star formation rate last- 
ing 10^ yr. More prolonged star formation episodes appear 
unlikely, given the very str ong feedback effects expected 
from pop HI starbursts (e.g. Ijohnson et "al]|2009l ). 

• Typical star formation efficiency (SFE), e. Here, 
the SFE is defined as the fraction of the baryons within the 
halo that is converted into stars throughout the duration of 
the pop III starburst (see above). In our model, the lumi- 
nosities of pop III galaxies simply scale with this parameter. 
For simplicity, we assume all halos to have a baryon fraction 
equal to the cosmic average f2bar/f^M, which results into the 
following relation between SFE and mass in pop III stars: 

i\i"popIII = eMhalo^^J^. (1) 

It should be noted, however, that many simulations pre- 
dict that halos with masses in the ~ lO'^-lO* Mq range 
may have their baryon fractions reduced by factors of a few 
due to radiative feedback f rom surrounding structures (e.g . 
iMesinger fc Diikstra ,20o"i: iRicotti. Gnedin fc Shull IbOOsi ) 



Mesmger &: Liiikstra ,2UUai: IKicotti. Unedm fe snull lll^UUBI I. 
In a recent paper, " iSafranek-Shrader et al.1 (|2012l ) present 



Yggdrasil model results available at: www.astro.su.se/~ez 
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simulations which suggest that the SFE of a pop III galaxy 
in a 3 X 10^ Mq halo at 2; ~ 12 may be logjQ e ~ —3. We 
therefore adopt this as our fiducial value, and investigate 
the circumstances under which pop III galaxies with this 
SFE would be detectable with the JWST. Given the sub- 
stantial uncertainties still attached to the typical SFE, we 
will also consider a wider range of values (up to logj^Q e = 0) 
when discussing the limits and detection prospects relevant 
for the HST. This covers the SF E range predicted in the 
simulations bv lWise fc CenI l|2009t ) (logj^Q e f» — 1 to —3) for 
a M ~ 10* M(T) halo and the de fault value (logi,, e ^ -2.2) 
assumed bv lPawlik et al.l l|201ll l for their pop III galaxies. 

• Lifetime of distinct pop III spectral signatures, 

Tpopiii. Even though the ejecta from the first pop III su- 
pernovae may take ~ 10^-10* yr to cool sufficiently to be 
used in the subse q uent formation of pop II and pop I stars 
(|Greif et al.l l2007l : iRitter et al.l 120121 '). the timescale during 
which pop III galaxies display unique spectral signatures 
may be much shorter. In principle, Tpopiii depends on the 
spectral signature considered, the IMF, /osc and the star 
formation history. For instanta neous-burst mode l s, and the 
colour signatures considered bv lZackrisson et al.l (|2011al lb[). 
Tpopiii is < 10^ yr, and may be as short as « 3 Myr if the IMF 
is extremely top-heavy (pop III.l) or even sho rter if feedback 
quick ly removes the gas from these systems jjohnson et al.l 
I2OO9I ). This parameter is used to reject pop III galaxies with 
ages in excess of Tpopiii from our final mock samples. An 
effective age truncation may, however, set in even sooner if 
the actual lifetime of a given pop III galaxy turns out to be 
lower than Tpopiii. In the case of Tpopiii = lO'^ yr, this hap- 
pens for our instantaneous-burst pop III.l galajcies, which 
remain luminous for just « 3 Myr before fading from sight. 

We adopt Tpopiii = lO'^ yr as our default value, but also 
consider TpopIII 3-8 short as 10^ yr. A more extended star 
formation history could potentially allow for a higher Tpopni, 
but given our maximum burst length of 10^ yr (see above), it 
is unlikely to be considerably longer than our default value. 
Moreover, the exact value of Tpopiii turns out to be one of the 
least critical parameters in our simulations (see Sect. 15.31) . 



All pop III fluxes are computed under the assumption 
of zero dust extinction. The flux at wavelengths shortward 
of - and including - Lya is moreover assumed to be zero to 
reflect the high level of absorption in the neutral intergalac- 
tic medium at the epochs considered {z ^ 7). Intrinsically, 
the LyQ emission line is expected to be one of the strongest 
features in the spectra of pop III galaxies with nebular 
emis sion, and could give rise t o distinct colour signatures 
|Peilo & Schaerer 2003; Zackri sson et al.l 12011^ ). While a 
certain fraction of the Lyo? photons could potentially be 
transmitted through the intergalactic medium (IGM) at 2; « 
7-8 d ue to outflows, source clustermg an d patchy reioniza- 
tion (|Diikstra. Mesinger fc Wvithd I2OIII : iDaval fc Ferraral 
I2OIII ). a high Lya transmission fraction would require very 
fortunate sightline and can most likely apply only to a small 
fraction of pop III galaxies in our mock samples. For this 
reason, we conservatively assume the flux of the Lya line to 
be zero, even though this may not always be the case. 



2.5 HST and JWST surveys 

The fluxes of pop III galaxies are predicted in two differ- 
ent fllter systems: the HST/WFC3 fllters used by the Clus- 
ter Lensing And Su pernova survey with Hubble {CLASH; 
IPostman etal] l2012l i and the JWST filters relevant for a 
hypothetical JWST survey of the J0717 cluster. 

CLASH is a 524-orbit multi-cycle treasury program tar- 
geting 25 galaxy clusters (including J0717) in a total of 16 
filters with the ACS and WFC3 instruments. The CLASH 
sample is larger and less biased than previous space-based 
cluster lensing surveys, with minimal lensing-based selec- 
tion that favors systems with overly dense cores. Specifi- 
cally, twenty CZ/j45'if clusters are solely X-ray selected to be 
massive {kT > 5 keV; 5 - 30 x lO" Mq) and, in most 
cases, dynamically relaxed. Five additional clusters were 
selected for their lensing strength (large Einstein radii or 
possibly high-magnification; to further quantify the lens- 
ing bias on concentration measurements, an d yield some of 
the h ighest-resolution dark matter maps; see lPostman et all 
|2012| for further details). For objects at 2 J5 7 (the range 
relevant for our simulations), detections are expected in at 
most 6 fihers (F850LP, F105W, FUOW, F125W, F140W 
and F160W), due to IGM absorption at wavelengths short- 
ward of Lya (A ^ 9700 A at these redshifts). The lOcr 
(5cr) detection limits in these filters are expected to be 
msso = 26.0 (26.7), mios = 26.6 (27.3), mno = 27.0 (27.8), 
mi25 = 26.5 (27. 2), mi4o = 26.7 (27.4 ), mieo = 26.7 (27.5) 
AB magnitudes l|Postman et al.|[2012l ). The field of view of 
HST/WFC3 (2.3' X 2.1') is too small to cover the entire 
J0717 field (Fig. [1] 4x4 arcmin), and when simulating the 
CLASH source counts we therefore just analyze an area of 
4.08 arcmin'^ towards the centre of the field, in rough agree- 
ment with how the CL^S-ff observations of this cluster have 
been carried out. This region is « 15% smaller than the full 
HST/WFC3 field of view, due to a rotation of the camera 
during the different exposures. 

Our hypothetical JWST multiband imaging survey is 
assumed to target J0717 with either tcxp = 3 h or tcxp = 30 h 
exposures per filter. The resulting IOct detection limits in all 
of the JWST/NIRCam and in the two bluest JWST/MIRI 
broadband filters are listed in Table [T] At similar wave- 
lengths (~ 1.50 pim), this survey reaches ~ 2.0 and 3.3 
mag (for toxp = 3 h and iexp = 30 h exposures respec- 
tively) fainter than the lOcr limits of CLASH. Of course, 
since JWST allows for observations at longer wavelengths 
than is possible with HST, the range of detectable object 
redshifts is also much larger. The field of view of NIRCam 
is f» 9.3 arcmin'^, divided into two 2.16 x 2.16 arcmin^ mod- 
ules with a ~ 50 arcsec gap in between. When simulating the 
NIRCam source counts in the J0717 field, we have therefore 
assumed a camera rotation that places the field of view di- 
agonally across the cluster (covering the lower left to upper 
right regions in Fig.[T}, thereby capturing most of the lensing 
caustics where the highest magnifications are expected. 

The pop HI galaxy number counts predicted in the 
J0717 field for the toxp = 3 h and 30 h detection thresh- 
olds are compared to what one would expect from ultra-deep 
field (UDF) style JWST observations of an unlensed field. 
In the latter case, the exposure times are assumed to be 
toxp = 100 h per filter, in rough agree ment with the deepes t 
-ffST observations carried out so far (jBeckwith et al.ll2006l ). 
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Figure 3. Predicted number of lensed pop III galaxy images as a function of JWST/F200W magnitude limit in the a) J0717 field and 
b) in an unlensed UDF. The coloured lines represent pop III galaxy models with the fiducial star formation efficiency (logj^Q e = —3) 
and different combinations of IMF and Lyman con tinuum e scape fraction: pop III.l, /esc = (red); pop III.l, /eac = 1 (orange); pop 
III.2, fesc = (green); pop III.2, /esc = 1 (cyan); iKroupal (jlooj), /esc = (blue). All models assume an instantaneous burst of star 
formation and Tpopm = 10^ yr. The dashed vertical lines indicate the Scr detection thresholds expected for different exposure times. 
Certain varieties of logj^Q e = —3 pop III galaxies are in principle detectable in the J0717 field (red, orange and green lines), whereas all 
such models lie well below the detection threshold of the UDF. 



The computational scheme used to simulate the JWST UDF 
is identical to the one outlined above, except that the pro- 
jection through J0717 magnification maps is omitted. This 
means that we implicitly assume the magnification across 
the UDF to be unity. One may argue that there is no such 
thing as a completely unlensed field at high redshift, since 
foreg round galaxies at z « 1-2 may magnify objects at 2 > 6 
ijWv ithc ot al. 20il|). For this reason, our results for the 
UDF should be considered somewhat conservative. Since we 
do not consider additional lensing boosts from intermediate- 
redshift galaxies located behind the J0717 cluster, a similar 
effect may be present in the J0717 field as well. 

In general, surveys of lensed fields can always reach 
deeper than blank-field surveys - and in just a fraction of 
the exposure time. However, blank-field surveys may still 
be more successful, depending on the luminosity function 
and the clus tering properties of the high-red shift objects 
targeted (e.g. lMaizv et al.ll2010l : lBouwens etHlbooa ) . High- 
redshift objects that are highly clustered or exhibit very low 
surface number densities may be impossible to detect behind 
a single lensing cluster, even if the magnification is sufficient 
to lift the intrinsic fiuxes of this class of objects above the 
survey detection threshold. Luckily, the surface number den- 
sities of pop HI galaxies is hig h (Fig. [Sjl and their clustering 
weak llStiavelh fc Trentill2010l ). 



3 DETECTING POP III GALAXIES WITH 

JWST 

In Fig. |3^, we plot the expected number of lensed pop HI 
galaxy images A'^popiii in a single JWST / NIRcam J0717 
field, as a function of detection threshold in the F200W 
filter (the most sensitive NIRCam passband, allowing de- 
tections of objects up to 2 ~ 15). The coloured lines repre- 



sent a set of different pop HI galaxy models, all with SFE 
logjQ e = —3; (see Sect l2.4l for a discussion). The dashed ver- 
tical lines indicate the 5a detection thresholds corresponding 
to exposure times of toxp ~ 3 and 30 h per pointing. Detec- 
tions of pop HI galaxies are deemed possible whenever a 
model track ventures brightward of one of these limits (to 
the left in Fig. ^ and simultaneously predicts > 1 image 
(logjQ Afpopiii(m ^ miimit) > 0). The corresponding results 
for the JWST UDF are presented in Fig. (SJd. Since an un- 
lensed field probes a much greater background volume than 
a high-magnification field like J0717, the image counts are 
steeper in the UDF and reach higher maximum values, al- 
beit at much fainter fiux levels. As seen in Fig. pop HI 
galaxies with logj^Q e = —3 are detectable through JWST 
observations of a strong-lensing cluster like J0717, as long 
as the IMF is top-heavy (pop III.l or pop HI. 2) and the 
escape of Lyman continuum photons into the IGM is minor 
(red and green lines in Fig. [Sjl. Galaxies with a pop III.l 
IMF (characteristic mass ~ 100 Mq) are the easiest to de- 
tect, despite having the shortest lifetimes (« 3 Myr for an 
instantaneous burst). With this IMF, even pop HI galaxies 
with logj^Q e = —3.5 are potentially detectable, provided that 
fesc is low. This represents a hard limit - if pop HI galaxies 
typically attain SFEs lower than this, they will inevitably 
be beyond the reach of JWST at 2 > 7, since ultra-deep ob- 
servations of unlensed fields are predicted to fare even worse 
(Fig. |3Jd). In fact, none of our logjQ e = —3 pop HI galaxy 
models are able to make it above the tcxp ~ 100 h detection 
threshold of the JWST UDF. When a similarly top-heavy 
IMF is adopted in realistic simulations of ^ 10* Mq halos, 
strong feedback causes substantial Lyman continuum leak- 
age into the IGM, thereby violating the /esc = assumption, 
even for SFEs as low as logjQ e ~ —3 to —4 jjohnson et al.l 
l2009t ). Even if fesc = 0, one of our pop III.l, logj^g e — —3 
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Table 1. JWST lOcr detection limits (AB magnitudes) for 3h, 30h 
and lOOh exposures per filter based on the JWST exposure time 
calculator, version 1.4. The corresponding 5(t limits lie 0.75 
mag faintward of these thresholds. 



Instrument 


Filter 


A {^lm) 




ni30h 


™100h 


NIRCam 


F070W 


0.70 


28.0 


29.4 


30.1 




F090W 


0.90 


28.5 


29.8 


30.4 




F115W 


1.15 


28.6 


29.9 


30.5 




F150W 


1.50 


28.7 


30.0 


30.6 




F200W 


2.00 


28.9 


30.2 


30.9 




F277W 


2.77 


28.9 


30.1 


30.8 




F356W 


3.56 


28.7 


30.0 


30.6 




F444W 


4.44 


28.2 


29.4 


30.0 


MIRI 


F560W 


5.60 


25.7 


27.0 


27.6 




F770W 


7.70 


25.1 


26.4 


27.0 



scenarios does, however, remain marginally detectable in the 
J0717 field (orange line in Fig.EJi). 

The latest simulations do admittedly not give much sup- 
port for an IMF a s extreme as our pop III.l scenario (e .g. 
Stacy et all |2010|: IClark et all l201ll : iGreif et all l2011al lbl: 
Hosokawa et al.l 201ll ). The pop III. 2 case (characteristic 
mass ~ 10 Mq) is in better agreement with current pre- 
dictions, and logj^Q e = —3 galaxies of this type should be 
detectable in the J0717 field provided that /osc is low (green 
line in Fig. [3^). Pop III galaxies wit h an IMF res embling 
that of galaxies in the local Universe (|Kroupall20dll ) remam 
undetectable at this SFE, regardless of the typical /esc (blue 
line in Fig.[3li). 

The mean redshift of pop III galaxies with non-zero 
JH^Sr/F200W fluxes in our mock samples is z « 11, but the 
mean redshifts of objects sufficiently bright to be detected 
above the J lyST detection limits are always lower than this. 
In general, mock samples based on parameters that generate 
intrinsically faint pop III galaxies (e.g. a Kroupa IMF, high 
/esc and low SFE) have lower mean redshifts than mock sam- 
ples with intrinsically bright objects (e.g. top-heavy IMF, 
low /esc and high SFE). The mean magnifications /i of all 
images with non-zero F200W fluxes in our NIRCam J0717 
field are /i ~ 10, but since high magnifications are required 
to make intrinsically faint objects detectable, the mean mag- 
nifications of images above the thresholds are also higher 
- especially so for samples with intrinsically faint pop III 
galaxies. For these reasons, the mean redshifts and magnifi- 
cations of logj^g e = —3 pop III galaxies that make it above 
the J W5r detection thresholds in the J0717 field have z « 9 
and p ^ 40. 



4 DETECTING POP III GALAXIES WITH HST 

Pop III gala xies at z ^ 7 with star for mation efficiency 
logio ~ -3 l|Safi-anek-Shrader et al.llioT^ ') are right at the 
expected detection threshold of JWST (Sect. [21), and well 
beyond the reach of current telescopes. However, if the SFE 
were significantly higher (logj^g > —2), objects of this type 
would in principle be detectable with the HST aX z ~ 7-10. 
In the following, we explore the prospects of detecting pop 
III galaxies in the ongoing HST/CLASH survey of lensing 



clusters, and derive the SFE limits imposed by existing HST 
data on unlensed fields. 

4.1 Detecting pop III galaxies in the HST/CLASH 
survey 

In Fig. 2^, we present the average number of pop III galaxy 
images A'popm predicted within the CLASH J0717 field as a 
function of detection threshold in the HST/WFC3 F160W 
filter. Formally, the CLASH survey reaches somewhat deeper 
in the FllOW filter (detection hmit mno = 27.0 at lOcr) 
than in F160W (mieo = 26.7 at IOct), but since objects start 
to appear as dropouts in the FllOW filter already at « > 8, 
our simulations still predict more pop HI galaxies in F160W. 
The pop HI galaxy models plotted (coloured lines) are the 
same as those discussed in the context of JWST detection 
(Fig. [3|, only shifted to a star formation efficiency an order 
of magnitude higher (to logj^g ^ ~ ~2). Most of these models 
lie just below the detection thresholds in the CLASH J0717 
field, and only pop HI galaxies with a pop III.l IMF and 
with negligible Lyman continuum leakage (red line) would 
be marginally detectable in these data. 

Since the J0717 cluster most likely exhibits the largest 
high-magnification area (~ 3.5 arcmin^ for > 10 at z > 6) 
among known cluster lenses, no other single cluster in the 
CLASH survey is likely to do any better. However, given 
the size of the CLASH survey, the prospects of detecting 
a pop HI galaxy in the full CLASH data set (25 clusters) 
are of course considerably better than when just the J0717 
images are used. To assess the maximum number of pop HI 
galaxies that may be detected in the full survey, we estimate 
the expected total area corresponding to high magnification 
(/i > 10) covered by the CLASH cluster sample. In doing 
so, we simply assume that the mass slope of the 20 relaxed. 
X-ray selected CLASH clusters will distribute around the 
mean value of dlogS/dlogS ~ —0.55, where E represents 
the surface mass density and 9 an angle in the sky. This 
is the typical dlogE/dlogS value seen for well-known re- 
laxed (lensing) clusters such as Abell 1703, Abell 383, Abell 
1689, C10024, among other examples. Since the majority of 
CLASH clusters is X-ray selected, so that their lensing power 
is not biased in that sense, we can assume that they would 
follow this (steeper) profile typical for relaxed clusters, which 
results in a more moderate magnification. These slope val- 
ues, according to models of several similar relaxed lensing 
clusters as mentioned above, yield a high-magnification area 
of ~ 1.55±0.4 arcmin^ per cluster on average, so that 20 such 
clusters should supply ~ 31 arcmin^ of area with /i > 10. 
For the five (expected) higher-magnification clusters we as- 
sume a shallower profile, as seen in our previous analysis of 
MACS0717 and MACS1149, or the full sample of 12 higher- 
z MACS clusters ijZitrin et al.ll2"oTll ). Average values therein 
are taken as representative for this part of the sample, with 
~ 2-2.5 arcmin'^ per cluster. In total, this yields a total area 
of ~ 42 arcmin^ for the CLASH sample. This is « 12 times 
greater than the corresponding area in J0717 alone, which 
means that the predicted pop HI galaxy number counts may 
be shifted upwards by approximately this factor once the 
full CLASH sample is considered. To illustrate the effects 
of this. Fig. 2)3 presents the image counts expected from 
the complete CLASH survey for the same logj^g e = —2 pop 
HI galaxy models as presented for just the J0717 field in 
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Figure 4. Predicted number of lensed pop III galaxy images as a function of HST/F160W magnitude limit in the a) CLASH J0717 
field and the b) full CLASH sample. The coloured lines represent pop HI galaxy models with star formation efficiency logjQ e = —2 and 
different combinations of IMF and Lyman continuum esca pe fraction: pop III.l, /osc = (red); pop HI. I, /esc = 1 (orange); pop HI. 2, 
/esc = (green); pop III. 2, /esc = 1 (cyan) : [Kroupal ll200fl ). /esc = (blue). All models assume an instantaneous burst of star formation 
and Tpopiii = 10^ yr. The dashed vertical lines indicate the 10a (thick dashed) and 5cr (thin dashed) detection thresholds of the CLASH 
survey. While no logj^Q e = —2 pop HI galaxies are expected in the J0717 data alone (a), a small number of logjg e = —2 pop III.l or pop 
HI. 2 galaxies with /esc = (red and green lines) may in principle turn up once the full CLASH data set is considered (b). 



Fig. [4^. Two of the pop III galaxy models (pop III.l and 
pop III. 2 with /esc = 0; red and green lines) are now shifted 
into the detectable range at logj^Q A'^popiii(m ^ mumit) ^ 0. 
Of course, the caveats concerning these parameter combina- 
tions discussed in Sect. [3] apply here as well. 

The mean redshift of all pop III galaxies with non-zero 
fluxes in the HST/F160W filter is z « 10, but in analogy 
with the case for JWST, the objects that actually appear 
above the detection thresholds tend do be at somewhat lower 
redshifts. For the models depicted in Fig. 3] the objects ap- 
pearing on the brightward side of the detection thresholds 
have z ~ 8-8.5. The mean magnifications p. of these images 
is /i « 40-60. 



4.2 Constraints set by the UV luminosity function 

While many combinations of model parameter values (see 
Sect. 12.4")) produce pop III galaxy models with fluxes far be- 
low the detection thresholds of both the JWST (Fig.O and 
HST (Fig. 3)), some mock samples with very high SFEs result 
in more detectable sources per surface area than observed 
in the HST/WFC3 data on the Hubble Ultra Deep Field 
(HUDF) and the adjacent Great Observatories Origins Deep 
Survey (GOODS) fields. The parameter combinations pro- 
ducing such violations can therefore be ruled out. The up- 
per limits on the typical SFE resulting from a comparison of 
the pop III luminosity functions (LFs) to th e observational 
constraints on the LF at z ~ 7, 8 and 10 (|Bouwens et al.l 
l201ll : IOesch et allbOlj ) are hsted in Table. [1 For the very 
brightest pop III galajcy models, the resulting upper limit is 
logj^Q e < —1.5, whereas for the very faintest models, the LF 
is unable to set any meaningful upper limits on the SFE. 

The basic principle is illustrated in Fig. [5] for an 
instantaneous-burst model with a pop III. 2 IMF, /esc ~ 



and Tpopiii = 10^ yr. A star formation efficiency of logjQ e = 
— 1 (red line) would overshoot the observed rest-frame ul- 
traviolet LF a.t z — 7 (left panel), and can therefore be 
ruled out, even though it is marginally consistent with the 
LF at z = 10 (right panel). An SFE of log^o e = -1-5 (or- 
ange lines) can on the other hand not be ruled out at any 
of these redshifts. In this case, a detailed calculation places 
the upper limit on the SFE at logj^g e ^ —1.2 (see Table 
[2|. This procedure does not in any way address the issue of 
whether pop III galaxies may already exist in current data 
(Zackrisson ct al. 2011b), but simply ensures that clearly 
unrealistic models are not considered any further. In fact, 
models located right at the SFE limits presented in Table. [2] 
would be consistent with having a few pop III galaxies in 
current HUDF and GOODS samples at z » 7-10. In most 
cases, the strongest constraints are set by the LF at 2: ~ 7, 
since this is where the observed LF reaches the deepest in 
terms of intrinsic luminosity. This technique of using ob- 
served source counts to constrain the pro perties of pop III 
galaxies is similar in spirit to that used by Zackrisso n et al.l 
(2010b) to impose limits on high-redshift supermassive dark 
stars. Such objects could potentially attain very high sur- 
face number densities, but in practice violate the observed 
LF unless the model parameters are tuned to make these 
objects exceedingly rare. 

It is worth noting that many of the models that remain 
viable once the LF limits have been applied still predict 
sufficiently many pop HI galaxies to dominate and alter the 
slope of the overall LF a few magnitudes below the current 
detection threshold. This is especially the case at z « 10, 
due to the large volume densities of pop III galaxies at this 
redshift. In Fig. [Bl we demonstrate this by plotting the z « 
10 LFs resulting from a number of different pop III galaxy 
models with fixed logj^g e — —1.5. Hence, if the typical pop 
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Figure 5. Example of how the observed rest-frame ultraviolet luminosity function (LF) at high rcdshifts constrains the typical star 
formation efficiencies of pop III galaxies. The panels illustrate how the observed LFs (Bouwens ct al. 2011; Ocsch et al. 2013) s.t z ^ 7 
(left panel) and 2 10 (right panel) set upper limit to the star formation efficiency e of pop III galaxies with a pop III. 2 IMF, /osc = 0, 
■'"poplll = 10^ yr and an instantaneous burst of star formation. The different line colours represent different e values. In this case, the 
logj^Q e = —1 model (red line) overshoots the faintest bin of the observed LF at 2 fii 7 (left panel) and is therefore ruled out, whereas 
models with logj^Q e = —1.5, —2.0, —2.5 and —3.0 (orange, green, cyan and blue lines) remain viable. A more detailed analysis places the 
upper limit at logjQ e < —1.2 (in between the models represented by the red and orange lines). Adopting a model with e close to this 
limit would, however, imply that a significant fraction of the faintest objects detected at 2 = 7-10 in fact are pop III galaxies. 



Ill SFE were this high, future extensions of the LF to fainter 
magnitudes would be expected to detect an abrupt upturn. 

Changing the IMF from a very top-heavy form (pop 
III.l) to a more moderatel y top-heavy v ersion (pop III. 2) 
or a perfectly normal IMF l|Kroupall200ll ) has the effect of 
increasing the lifetimes of the model galaxies, at the expense 
of the lowering their peak luminosities, thereby shifting the 
predicted pop III galaxy luminosity distribution upwards 
and to the right in these diagrams (compare red and orange 
lines in Fig. [S] or the orange line in Fig. [S] to the blue line 
in Fig. [6|. Adopting a more extended star formation history 
than an instantaneous burst has a very similar effect, since 
this also lowers the peak luminosity but increases the life- 
time (compare the orange line in Fig. [5] to the cyan line in 
Fig. [S]). Since pop III galaxies typically attain their greatest 
luminosities at very low ages, adopting a Tpopm lower than 
lO'^ yr primarily results in a truncation of the low- luminosity 
tail of LF, with a pronounced effect on the overall shape of 
the pop III LF but with little relevance for current LF mea- 
surements. 

For a number of reasons, the star formation efficiency 
limits presented in Table [2] should be considered approx- 
imate only. Whereas the absolute magnitudes plotted in 
Fig. [5] and [6] are based on the same HST/WFC3 filters as 
the observations (F125W for « « 7; F160W for z f» 8 and 
z « 10), the mock sample data are based on simulations 
statistics within redshift bins of depth Aa = 1.0 centered at 
2: = 7, 8 and 10, with all pop III galaxy magnitudes com- 
puted using the central redshift. This is admittedly a crude 
approximation, since the unusual SEDs of pop III galaxies 
may result in a redshift selection function that is very dif- 
ferent from that of more mundane galaxies. This is particu- 
larly the case if /Lva is non-zero for some of these galaxies 
l|Zackrisson et al1l2011bl ). Also, the results presented in Ta- 
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Figure 6. Examples of how the predicted pop III galaxy lumi- 
nosity function (LF) at 2 10 is shaped by the pop III galaxy 
model parameters. The differently coloured lines represent the 
LFs predicted in the case of a pop III.l IMF with /osc = 1 (red 
line); a pop III. 2 IMF with /esc = 1 (orange); a pop III. 2 IMF 
with /osc = and Tpopin = lO'' yr (green); a pop III. 2 IMF with 
/osc = and a lO'' yr burst (cyan); a pop III. lKroupal 1I2OOII ) IMF 
with /osc = (blue). The black symbols represent the observed 
rest- f rame ultraviolet LF constraints at this redshift llOesch et al.l 
All models assume logj^Q e = — 1.5 and remain viable once 
the LF constraints (Table [2]l have been applied. For this star for- 
mation efficiency, the overall LF is expected to undergo dramatic 
changes a few magnitudes below the current detection threshold 
due to the large volume densities of pop III galaxies at this red- 
shift. 
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Table 2. Constraints on the star formation efficiency e of pop 
III galaxies derived from the observed rest-frame UV luminosity 
functions at 2 = 7-10 
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ble [2] do not take the LF contribution from more mundane 
galaxies (which are Ukely to be far more abundant than pop 
III galaxies) into account, which makes the upper limits on 
the star formation efficiencies conservative. 



5 DISCUSSION 

5.1 Lensed versus unlensed fields 

In Fig. [3 we compare the expected performance of both 
HST/CLASH and our hypothetical JWST J0717 survey to 
that of ultradeep (texp = 100 h) JWST/NIRCam observa- 
tions of an unlensed field, by plotting the total number of 
pop III galaxy images A'^popiii above the survey 5<j detec- 
tion thresholds as a function of the pop III galaxy SFE. 
The results in this figure are based on a model with a Pop 
III. 2 IMF, /esc = and Tpopiii = 10^ yr, but the relative 
difference between the survey results are fairly similar for 
all models considered. Since lensed fields cover a smaller 
background volume, more pop III galaxies are predicted in 
the JWST UDF at high SFEs, but the numbers drop more 
rapidly with decreasing SFE than in the JW^ST J0717 obser- 
vations. For this particular pop III galaxy model, only JWST 
observations of J0717 are able to reach logj^g ^ ^ —2.5 with 
A'popiii > 1 (i.e. logj^g A'^popiii > 0). In relative numbers, a 
JWST survey of J0717 probes SFEs 0.4 dex (3 h exposures) 
and 1.0 dex (30 h exposures) deeper than the JWST UDF. 

Included in Fig. [7] are also the pop III image counts 
predicted for the HST/CLASH survey. A direct compari- 
son between the HST and JWST results for J0717 indicates 
that the JVFST survey of this cluster would be able to probe 
SFEs (or, equivalently, intrinsic luminosities) that are lower 
by ~ 1.2 dex (3 h exposures) and ~ 1.8 dex (30 h expo- 
sures). When compared to the expected performance of the 
fuU HST/CLASH sample, this JiyST survey instead reaches 
^ 0.7 dex (3 h exposures) and « 1.3 dex (30 h exposures) 
deeper in terms of SFE. In summary, observations of a high- 
magnification cluster like J0717 should offer the best chances 
of detecting pop III galaxies with JWST. 
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Figure 7. The number of pop III galaxy images predicted above 
the 5cr detection thresholds in different surveys, as a function of 
the pop III star formation efficiency e. The blue lines represent 
the expected results from the HST/CLASH images of the J0717 
galaxy cluster (blue solid) and of all the CLASH clusters (blue 
dashed) in the F160W filter. The red lines indicate the corre- 
sponding results from JWST F200W images with 3h (red solid) 
and 30h (red dashed) exposure times of J0717. These lensing sur- 
veys are contrasted with the corresponding number of pop III 
galaxies predicted for a JWST/UDF lOOh survey (black line) of 
an unlensed field. A model with a Pop III. 2 IMF, fesc = and 
Tpoplll = 107 yr is assumed. The gray region represents the SFE 
range ruled out by existing LF constraints fSect. |4?2] |. While the 
JWST/UDF images are expected to produce far more pop III 
galaxies than the J0717 observations at high star formation ef- 
ficiencies, only J0717 observations can probe pop III galaxies at 
logio ^ -2-5. 

5.2 Image multiplicity 

In situations where our lensing simulations predict 
logiQ Afpopiii(m ^ miimit) > (i.e. more than one detectable 
image), the number of lensed pop III galaxy images will 
typically not equal the number of unique pop III galaxies 
detectable in the surveys we consider, since a fraction of the 
images will stem from the same source. We find that the 
image multiplicities of our mock samples usually lie in the 
range « 1.2-2. Hence, our A'popiii estimates in Figs. |31 HI 
and [7] should be divided by a factor in this range to recover 
the number of unique pop III galaxies. In the case of the un- 
lensed JWST/UDF Held (Fig.lSJa), no correction is of course 
required. The reason why we present the number of images, 
rather than the number of unique galaxies, is that having 
more images tends to be advantageous - some of the im- 
ages (both in lensed and unlensed surveys) may be rendered 
useless due to blending with foreground objects (Milky Way 
stars and low-redshift galaxies). 

In principle, the image multiplicity depends on the po- 
sitions of the sources with respect to the caustics of the 
lensing cluster. However, our simulations also indicate a cor- 
relation between image multiplicity and the magnifications 
of the images that make it above the detection threshold. 
This makes the typical image multiplicities sensitive to the 
intrinsic luminosities of the pop HI galaxies which, in turn, 
depend on the model parameters described in Sect. 12.41 The 
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fact that the image multiphcity in many cases is < 2 indi- 
cates that not all lensed pop III galaxies produce detectable 
multiple images. There are several reasons for this. To begin 
with, not all highly magnified pop III galaxies are multiple 
imaged - our lens mapping predicts a significant fraction of 
single images close to but outside the caustics (e.g. « 50% of 
all images with ^ > 20). Secondly, the flux ratios of the mul- 
tiple images are sometimes too large to allow for more than 
one to make it above the detection threshold - flux ratios of 
up to a factor of ~ 2 (corresponding to ~ 0.75 magnitudes) 
between the brightest and second brightest images are not 
uncommon. Finally, some multiply imaged systems display 
such large image separations (> 1') that even though several 
images are in principle detectable in the J07I7 field (Fig.[T]), 
only a single one appears within the part of the cluster that 
we assumed to be covered by the observations. 

5.3 Star formation history and duration of the 
pop III phase 

The results presented throughout Sect. [3] and 14. II are based 
on the assumption that pop III galaxies form their stars 
more or less instantaneously and regain their pop III spectral 
signatures for Tpopiii = 10^ yr. 

Contrary, perhaps, to naive expectation, the Tpopiii pa- 
rameter is not necessarily of critical importance for the 
prospects of detecting pop III galaxies in the surveys we 
consider. While reducing rpopiii from 10'' yr to 10*^ yr for- 
mally reduces the number densities of all pop III galaxies by 
an order of mag nitudeS at any redshift, the numbers of ob- 
jects that appear above the relevant detection thresholds in 
figures like Fig.|3]and |4]typically change by less than « 50%. 
This happens because it is only the very brightest pop III 
galaxies that have any sporting chance of being detected. 
Since pop III galaxies experiencing an instantaneous burst 
of star formation attain their peak luminosities already in 
K. 1-2 X l(f yr (depending on exact IMF and filter), remov- 
ing objects with ages above this limit has little or no impact 
on the high-luminosity tail of the intrinsic pop III luminos- 
ity distribution. This effect can be readily seen by compar- 
ing the intrinsic luminosity functions for pop III. 2, /esc = 0, 
logj^Q e = —1.5 galaxies with Tpopiii = lO'^ yr (orange line in 
Fig. O and Tpopiii = 10*' yr (green line in Fig. [SJ - the two 
distributions are similar at the high-luminosity end but very 
different at lower luminosities. Gravitational lensing can in 
principle overthrow this trend by assigning high magnifica- 
tions to intrinsically faint galaxies and low magnifications 
to intrinsically bright galaxies, but in our simulations, this 
mechanism is not sufficiently efficient to make Tpopiii a cru- 
cial parameter in the current context. 

Extending the star formation history has the effect of 
lowering the peak luminosity per unit star formation effi- 
ciency and to prolong the brightest phase in the life of these 
objects. This makes all simulated pop III galaxies less lumi- 
nous but more numerous (as long as Tpopiii is not shorter 
than the duration of the starburst). When it comes to the 
prospects of detecting rare objects in lensed fields, boost- 
ing the intrinsic (unlensed) object number densities could 

* For pop III.l galaxies, which have total lifetimes of as little as 
K, 3-4 Myr, the correction factor is closer to 1/3 



in principle be advantageous, but in our simulations it does 
more harm than good, partly because we impose a ^ 100 
threshold. When going from an instanstaneous burst to a 
lO' yr star formation episode, we find that the SFE on av- 
erage needs to be boosted by Alogj^ge ~ 0.2-0.4 dex (i.e. 
a factor of « 2-3) to produce the same number of objects 
above the detection thresholds in our lensed HST and JWST 
surveys. 

5.4 Uncertainties in the lensing maps 

The IZitrin et al] l|2009l 'l J0717 magnification maps used in 
our simulations come with a certain level of uncertainty, and 
this translates into an uncertainty in our estimates of the 
number of detectable pop III galaxies in the HST and JWST 
surveys of the cluster we consider. Since the deflection fleld 
scales with the lens-source (-Dis) to source (-Ds) angular-size 
distance ratios, cluster lensing maps require knowledge of 
at least one multiply-lensed source redshift to calibrate the 
mass model, and of other redshifts in order to deduce the 
mass profile and accurately assess the magnification. How- 
ever, by the time our J0717 magnification maps were pro- 
duced, no such spectroscopic redshift s were available. Ba sed 
on only two available imaging bands, IZitrin et al.l l|2009l ) in- 
stead estimated a redshift of 2; « 2-2.5 for the main source to 
which their model was calibrated (system #1 in their nota- 
tion), whereas other dropout multiply-lensed objects found 
in their work were used to pin down the pro file and the re- 
sulting magnification. The lZitrin et al.l l|2009t ) lensing maps, 
and by extension our simulation results, are therefore im- 
plicitly based on the assumption that the redshift of system 
#1 is z « 2. Recently, the first spectroscopic information 
was published for objects behind clu ster, revealing that # 1 
in fact is at a redshift oi z — 2.963 l|Limousin et al.ll201ll ). 
This entails two possible systematic errors which we now 
make an effort to assess. 

First, the lensing distance ratios for higher-z source 
planes should now be slightly lower. Luckily, due to the con- 
verging nature of the Di^/Ds for higher source redshift for 
a fixed lens, the corrections for the higher-redshift source 
plane angular-size distances are minor and decrease with in- 
creasing source redshifts. For example, assuming that #1 
is at z = 2 as we did, results in a relative Di^/Ds factor 
of ~ 1.3 for a source-plane redshift of Zg = 10. By instead 
using the correct redshift of Zs — 2.963, the relative ratio 
becomes « 1.2. To assess the impact of this on our results, 
we have produced magnification maps for a typical redshift 
of Zs = 10 according to both #1 redshift options and com- 
pared the magnification properties. We find that the area of 
/I > 10 we used in our calculations is only 8% larger; and for 
a threshold of ^ > 20 or higher, larger by < 1%, therefore 
rendering a possible systematic bias negligible. 

Second, due to th e higher redshifts of system #1 and #5 
than was assumed in IZitrin et al.l l|2009l ) , the cluster mass 
profile might also be affected at some level. Due to the ex- 
tended merging nature of this cluster, the profile is still ex- 
pected to be very shallow, so that any impact on the mass 
profile is correspondingly expected to be very small. To ex- 
amine this effect, we have checked that within the given 
uncertainties, our model is commensura te with the new red - 
shift measures for the systems found in lZitrin et al.l l|2009l ). 
especially systems #3 and #5 which expand the redshift 
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range, and reproduces these systems at their given redshifts. 
We find a general good agreement, although we note that 
the mass profile of the model in use may be slightly shallower 
in the very inner core than needed to perfectly match the 
Dis/Ds growth. We therefore compare the high magnifica- 
tion areas with a variety of steeper models, and find that the 
size of the high-magnification areas are subject to a < 10% 
uncertainty because of this. 

Since the expected number of detectable pop III galax- 
ies is proportional to the size of these high-magnification 
regions, we conclude that the n ew redshift information pro- 
vided by iLimousin et al.1 (|201ll ) would result in very small 
adjustments of our predictions. A further refinement of the 
mass model will soon also be enabled by the CLASH 16-band 
imaging of J0717 and the resulting photometric redshifts for 
its multiply-lensed systems. 



5.5 The colour signatures of pop III galaxies 

Throughout this paper, we have explored the pop III galaxy 
properties required to place them above the fiux detection 
thresholds of different HST and JWST surveys. However, 
having a galaxy bright enough for detection does not neces- 
sarily mean that one will also be able to identify it as a pop 
III object based on its colours. In fact, only certain subtypes 
of pop III galaxies are likely to display colour signatures that 
set them apart from from more mundane objects. 

Like all galaxies in the reionization epoch (z > 6), pop 
III objects will appear in dropout-selected samples because 
of absorption shortward of the Lya line by the largely neu- 
tral IGM. However, additional colour criteria are also re- 
quired, and these depend on the exact properties of the 
pop III galaxies. Pop III galaxies dominated by direct star 
light can be uniquely identified for ages at least up to « 10 
Myr in HS T or Jiygrfihers because of their very blue UV- 
slope s (e.g.lSchaerer fc Pelldll2005l : iRaiter et al]|2010l : llnoud 



I2OIII : IZackrisson et al.] 2011a). Pop HI galaxies dominated 
by nebular emission (all our /eac = 0, Tpopiii ^ 10^ yr mod- 
els) can be identified from their JWST colours because of 
their very stron g hydrogen Balmer emission lines bu t lack of 
[OIII] emission (lnouell201ll : IZackrisson et al.ll2011al '). Unfor- 
tunately, the colour signatures arising from this are unique 
only over very limited redshift intervals. The [OIII] line at 
5007 A also redshifts out of the wavelength range of the 
JWST/NIRCam instrument at 2 > 8, and even though the 
JWST/MIRI camera, can be used to probe this line at higher 
redshift, the lower sensitivity of this instrument makes this 
endeavour extremely ex pensive in terms of exposure time 
ijZackrisson et al1l2011al ). In general, we find that M/7?/ de- 
tections of objects appearing at the NIRCam detection limit 
(Fig. [31) would require such long exposure times (> 100 h per 
filter) that this approach scorns unrealistic. Pop III galax- 
ies with nebul ar emissio n can also exhibit extremely strong 
Lya emission l|Schaereill2002l . l2003l : iRaiter et al.il2blol '). and 
this can in principle produce distin ct colour signatures in 
either HST jZackrisson et al.ll2011bl ) or JWST filters. How- 
ever, since Lya is readily absorbed in the increasingly neu- 
tral IGM at 2; > 6, this requires a fortunate sightline and can 
most likely apply only to the objects at the lowest redshifts 
considered in our simulations. 

Realistically, colour criteria can probably only be used 
to select pop III galaxy candidates, whereas follow-up spec- 



troscopy will be required to achieve a more robust iden- 
tification of pop III objects, based either on UV contin- 
uum slope or the strength of emission lines like Lya, Hell, 
H/3, H a and [QUI] (iSchaerer 2002 2003: Schaerer fc Pellol 
200 d: IRaiter et al.1 120101 : llnouej I2OI1I : ICai et al.1 120111 
Zac krisson et al.l 1201 Ibl ) Fortunately, objects appearing 
above the detection thresholds of HST/CLASH (mAB < 27) 
are within range of existing near-IR spectrometers on 8- 
10 m groundbased telescopes (e.g. ISchenker et al.l |2012| : 
iLehnert et"ai]|2010l '). and objects appearing above the detec- 
tion limits of a lensed JWST survey (mAB ^ 29 for texp = 3 
h) are within the range of the JWST/NIRSpec spectrometer 
or planned groundbased telescopes like the Giant Magellan 
Telescop^ the Thirty Meter Telescop^ and the European 
Extremely Large Telescop^ The same exercise would be far 
more difficult for pop III galaxy candidates with compara- 
ble SFEs detected in an unlensed field like the JWST/UDF. 
For example, the lowest SFE at which the JWST/UDF still 
produces competitive pop HI galaxy counts compared to a 
lensed JWST survey (see Fig. [7|) converts into an unlensed 
pop HI galaxy fiux of rriAB ~ 30.5-31.0, which would be 
very challenging to probe spectroscopically in the foresee- 
able future. 



5.6 Hybrid pop III galaxies 

Our simulations trace the formation of metal-free halos that 
are likely to host pure pop III galaxies, but due to incom- 
plete mixing of metals in the interstellar medium, pop HI 
star formation may plausibly continue at some minor level 
within obj ects that have already started to form pop I/II 
stars (e.g. I Johnson et al.l l2009l : Isalvaterra et al ][20l3). In 
most cases, this will dilute any pop III spectral or colour 
signatures beyond recognition, but certain diagnostic crite- 
ria, like the strength of the Lya line, can in principle ap- 
ply even if just a minor fraction of the overall stellar mass 
of a high-redsh ift galaxy is in the form of pop III stars 
(|Zackrisson et a l. 2011b). Such hybrid pop HI galaxies are 
not included in our source count estimates, and could boost 
the prospects of detecting pop HI galaxies in both lensed 
and unlensed fields. 



6 SUMMARY 

Using simulations of the fiuxes of pop III galaxies, we argue 
that the observed galaxy luminosity function at 2: = 7 — 10 
can be used to place upper limits on the typical star for- 
mation efficiencies of these objects, and present such limits 
for a wide range of model parameters. By projecting these 
simulations through the magnification map of the galaxy 
cluster J0717 (arguably the best lensing cluster currently 
known for high-redshift surveys), we predict that the on- 
going HST imaging survey CLASH, which targets a total 
of 25 galaxy clusters - including the J0717 - potentially 
could turn up a small number of pop III galaxies if the typ- 
ical star formation efficiency is ~ 10~^. We also explore the 



^ http://www.gmto.org/ 

^ http://www.tmt.org/ 

^ http:/ /www. eso.org/sci/facilities/eelt/ 
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prospects of surveying the J0717 cluster with the upcoming 
JWST telescope, and find that this strategy would be able 
to detect pop III galaxies with star formation efficiencies an 
order of magnitude lower (~ fO~^ of the baryons converted 
into pop III stars). The gravitationally enhanced fluxes of 
such objects also make them potentially suitable targets for 
follow-up spectroscopy, whereas this would be exceedingly 
difficult for corresponding objects in unlensed fields. 
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